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Abstract. We construct a semiclassical Friedmann-Lemaitre- Robertson- Walker (FLRW) 
cosmological model assuming a running cosmological constant (CC). It turns out that 
the CC becomes variable at arbitrarily low energies due to the remnant quantum effects 
of the heaviest particles, e.g. the Planck scale physics. These effects are universal in the 
sense that they lead to a low-energy structure common to a large class of high-energy 
theories. Remarkably, the uncertainty concerning the unknown high-energy dynamics 
is accumulated into a single parameter v, such that the model has an essential predic- 
tive power. Future Type la supernovae experiments (like SNAP) can verify whether 
this framework is correct. For the flat FLRW case and a moderate value v ~ 10 -2 , we 
predict an increase of 10 — 20% in the value of $7 a at redshifts z = 1 — 1.5 perfectly 
reachable by SNAP. 
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Introduction 

The number and versatility of publications concerning the Cosmological Constant (CC) prob- 
lem, show that it is a fascinating interface between Cosmology and Quantum Field Theory (QFT). 
One can identify three distinct aspects of this problem: I) The famous "old" CC problem (see pQ) 
is why the induced and vacuum counterparts of the CC cancel each other with such a huge preci- 
sion; II) The coincidence problem (see the reviews pj) is why the observed CC in the present-day 
Universe [Hj is so close to the matter density p ; III) The dynamics responsible for the CC prob- 
lem: is the CC some form of dark energy associated to a special new entity? If future astronomical 
observations (like the SNAP project [1]) would discover that the CC depends on the redshift pa- 
rameter z, one could wonder whether this variation can be achieved only through the mysterious 
quintessence scalar field searched by everyone. The idea of a scalar field that may adjust itself 
such that the cosmological constant is zero or very small, has been entertained many times in the 
literature and since long ago jSJ E] . But in practice one ends up with all sorts of (more or less 
obvious) fine-tunings. As emphasized in 003111]) the quintessence proposal for problem III does 
not help much in solving problems I and II. In our opinion, before introducing an ad hoc field we 
had better check whether the variation of the CC can be attributed to a simpler, more economical 
and truly robust QFT concept, like e.g. the renormalization group (RG). 

The relation between the RG and the CC problem had been in air for a long time lllj . 
but only recently it has been settled into the phenomenological framework jSJ 1121 113j . A con- 
sistent formulation of the approach within the context of particle physics in curved space-time 
has been presented in [5], where the relevant RG scale fi has been identified with the energy 
of the gravitational quanta; that means, in the cosmological arena, with the Hubble parameter: 
fi = H. Other possibilities have been considered in the literature and in different contexts 3 . At 
first sight, the tiny present-day value Hq ~ 1.5 x 10 -42 GeV makes the RG running of the CC 
senseless, because all the massive fields should decouple. Of course, this is true if we apply the 
"sharp cut-off" scheme of decoupling 8 j. However, in ^2] a pertinent observation concerning the 
decoupling has been made. In fact, the decoupling is not sharp, and the quantum effects of a 
particle with mass m at the energy \x <C m, even if suppressed by the factor (^/m) 2 , can be im- 
portant in the CC context. Recent investigations of the decoupling in the gravitational theory |14j 
have confirmed explicitly this decoupling behavior for the higher-derivative sector of the vacuum 
action. However, the running of the CC and the inverse Newton constant is not visible within the 
perturbative approach of 14 , while alternative covariant methods of calculations are incompatible 
with the mass-dependent renormalization scheme and their development may require a long time. 
In this situation it is reasonable to apply the phenomenological approach. In the present paper we 
suppose that the decoupling of the CC occurs in the same way as in those sectors of the vacuum 
action, where we are able to obtain this decoupling explicitly [Tl|. As we are going to show, this 
approach may shed light to the coincidence problem II and leads to very definite predictions for 

3 In Ref.[8] the alternative possibility associating fi with p 1 ^ 4 is put forward, and in 1121 it is amply exploited. In 
Ref.|13) the RG scale /x was identified with the inverse of the age of the universe at any given cosmological time, i.e. 
/i ~ 1/t. This is essentially equivalent to our choice 0, because H ~ 1/t in the FLRW cosmological setting. 
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the variable CC, which can be tested by the rich program of astronomical observations scheduled 
for the next few years. 

Renormalization Group and Planck Scale Physics 

Consider a free field of spin J, mass Mj and multiplicities (n c , nj) in an external gravitational 
field - e.g. {n c ,ni/ 2 ) = (3,2) for quarks, (1,2) for leptons and (n c , no,i) = (1,1) for scalar and 
vector fields. At high energy scale, the corresponding contribution to the (3- function dA/dlnfi for 
the CC is 

(-1) 2J 

/3 A ( M >Mj) = ^—l— (J +1/2) n c nj Mj. (1) 

At low energies (/i <C Mj) this contribution is suppressed due to the decoupling. At an arbitrary 
scale fj, (whether greater or smaller than Mj ) , the contribution of a particle with mass Mj 
should be multiplied by a form factor F([j,/Mj) . At high energies F(fi 3> Mj) ~ 1 because there 
must be correspondence between the minimal subtraction scheme and the physical mass-dependent 
schemes of renormalization at high energies. In the low-energy regime \i <C Mj one can expand 
the function F into powers of njMj: 




Two relevant observations are in order. First, the term n = must be absent, because it would 
lead to the non-decoupling of Mj, with untenable phenomenological implications on the CC value. 
Indeed, for those terms in the vacuum action where the derivation of the function F(/i/Mj) is 
possible [Tl], the n = terms are really absent. Second, due to the covariance the number of 
metric derivatives (resulting into powers of H) must be even, and so there are no terms with 
odd powers in the expansion (J2J) - remember we identify fi with the Hubble parameter. No other 
restrictions for the coefficients k n can be seen. Indeed, in the H <C Mj regime the most relevant 
coefficient is k\ . In the rest of this article we develop a cosmological model based on the hypothesis 
that this coefficient is different from zero. This phenomenological input does not contradict any 
known principle or law of physics and at the present state of knowledge its validity can be checked 
only by comparison with the experimental data. 

As far as we suppose the usual form of decoupling for the CC, all the contributions in Q are 
suppressed by the factor of (fi/Mj) 2 and we arrive at the expression 

(-1) 2J 

/3 A Ou « Mj) ~ h y -j^T (J + 1/2) n c nj Mj fi 2 . (3) 

At the very low (from the particle physics point of view) energies [i = Hq ~ 1.5 x 10"~ 42 GeV , 
the relation (/i/Mj) 2 <C 1 is satisfied for all massive particles: starting from the lightest neutrino, 
whose presumed mass is m u ~ 10 30 i^o, up to the unknown heaviest particle M + < Mp. Then, 
according to ©, the total /^-function for the CC in the present-day Universe is, in very good 
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approximation, dominated by the heaviest masses: 

M+ 

01= Y, Pa(Mj) * j—rt oM 2 f. (4) 

Mj=m v ^ ^ 

Here we have introduced the following parameters: M is a mass parameter which represents the 
main feature of the total /3-function for the CC at the present cosmic scale, and a = ±1 indicates 
the sign of the CC /3-function, depending on whether the fermions (a = —1) or bosons {a = +1) 
dominate at the highest energies. The quadratic dependence on the masses Mj makes (3\ highly 
sensitive to the particle spectrum near the Planck scale while the spectrum at lower energies has 
no impact whatsoever on (3\. 

Having no experimental data about the highest energies, the numerical choice of a M 2 is 
model-dependent. For example, the fermion and boson contributions in © might cancel due to 
supersymmetry (SUSY) and the total /^-function becomes non-zero at lower energies due to SUSY 
breaking. In this case, the value of M 2 depends on the scale of this breaking, and the sign a 
depends on the way SUSY is broken. In particular, the SUSY breaking near the Fermi scale leads 
to a negligible j3\, while the SUSY breaking at a GUT scale (particularly at a scale near the Planck 
mass) provides a significant (3^. 

Another option is to suppose some kind of string transition into QFT at the Planck scale. Then 
the heaviest particles would have the masses comparable to the Planck mass Mp and represent 
the remnants, e.g., of the massive modes of a superstring. Of course, this does not contradict 
SUSY at the lower energies, when these string WIMPzillas decouple from the matter sector. But, 
according to Eq. (|4"|). they never decouple completely from the RG in the CC case. Below, we shall 
work with this option and take, for simplicity, M 2 = Mp. Let us clarify that this choice does not 
necessary mean that the relevant high energy particles have the Planck mass. The mass of each 
particle may be indeed smaller than Mp, and the equality, or even the effective value M > Mp, 
can be achieved due to the multiplicities of these particles. With these considerations in mind, our 
very first observation is that the natural value of the /3-function (JIJ at the present time is 

0i\ = j±y 2 M 2 ■ J = ^M 2 • H 2 ~ l^GeV\ (5) 

where c is some coefficient. For c = 0(1—10) the 0\ function is very much close to the experimental 
data on the CC [3]. This is highly remarkable, because two vastly different and (in principle) totally 
unrelated scales are involved to realize this "coincidence": Hq (the value of /i at present) and Mp, 
being these scales separated by more than 60 orders of magnitude! In this framework the energy 
scale of the CC is of the order of the geometrical mean of these two widely different ones, and 
only near Mp all scales become of the same order. Moreover at low energy the running of the 
CC proceeds smoothly and in the right ballpark required by the natural solution of problem II. In 
fact, this is not a complete solution of this problem. Of course Eq. Q holds at any cosmic scale, 
but the Friedmann equation leads to the relation H ~ \^A/Mp only at lower energies, while at 
higher energies the matter or radiation density dominates. Anyway, the coincidence between the 
supposed RG decoupling Q) and the Friedmann equation for the modern, CC-dominated, Universe 
looks rather intriguing and is worth exploring. 
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A cosmological FLRW model with running A 



Consider the implications for FLRW cosmological models coupled to Eq. (J3J). The first step 
is to derive the CC dependence from the redshift parameter z, defined as 1 + z = do /a, where ao 
is the present-day scale factor. Using the identification of the RG scale p, with H, we reach the 
equation 

dA 1 dH . oM 2 dH 

Tz ~ H ~dz~^ ~ ~dz" ( ) 

In order to construct the cosmological model, we shall use, along with Eq. Q, the Friedmann 
equation 

H 2 = - = — P + A - -, 7 
\ a J 3 a z 

where p is the matter/CDM/radiation density and the last curvature-dependent term can be 
presented as —k/a 2 = Hq£1° k (1 + z) 2 , where £l Q K is the spatial curvature parameter at present 
(z = 0), and can be written in terms of the usual cosmological parameters for matter and CC 
at the present time: = 1 — — il^. Also, the energy conservation law provides the third 
necessary equation 

where p is the matter /radiation pressure. One could wonder whether the non-local effects behind 
the renormalization group are compatible with the standard energy conservation equation Q . Let 
us remember, however, that the covariant form of the conservation law < V^T^ >= just reflects 
the covariance of the effective action and therefore does not depend on the non-localities which 
are always present in the quantum corrections. Furthermore, in a situation where the energy scale 
associated to the metric derivatives is very small, the leading effect of the CC scale dependence 
may be, according to our model, presented in a compact form due to the renormalization group. 
Then the proper form of the conservation law is the one we use here. 

As we shall consider both MD (matter dominated) and RD (radiation dominated) regimes, it 
is useful to solve the equations ©, Q, © using an arbitrary equation of state p = ap, with a = 
for MD and a = 1/3 for RD. The time derivative in © can be easily traded for a derivative in z via 
d/dt = — H (1 + z)d/dz. Hence we arrive at a coupled system of ordinary differential equations 
in the z variable. The solution for the matter-radiation energy density and CC is completely 
analytical. It reads a follows: 



and 



with 



p(z; v) = [p + — p° c ) (1 + zf - — pi (1 + zf (9) 
A(z;u)=A + ^-M 2 P (H 2 (z;u)-H 2 ) =A + p f(z) +p° c g{z), (10) 
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k(1 + 3a) , s (1 + z) C - 1 

g(z) = \ > z ( z + 2)+VK-± r4- r- (12) 

yy ' 2(C-2) v ; (l-^)(C-2) v ; 

Here po, p®, Aq and Hq are respectively the matter-radiation energy density, critical density, CC 
and Hubble parameter at present, and we have used the following dimensionless coefficients: 

a M 2 



12ttM| ' 



C = 3(1-1/) (a + 1), K=-2un° K . (13) 



It should be clear that there is only one single independent parameter in the model: v. In the 
limit v — > we recover the standard result for p{z) with constant CC (see, e.g. |15|). In order to 
avoid confusion, let us note that the above solutions for p{z\ v) and A(z; v) have no singularity in 
the limits £ = 2 and v = 1. It is easy to see that the function g{z) differs from zero only for the 
cases of non-vanishing spatial curvature k / 0, while another function f{z) remains significant 
even for the flat space case k = 0. Furthermore, in the limit v — ► , A(z; v) just becomes the 
z-independent standard CC. 

The formulas above represent the universal solution at low energies, when all massive particles 
decouple according to ©• In order to understand the limits of applicability of these solutions, 
let us notice that when the temperature of the radiation energy density achieves the Fermi scale 
T ~ Mp ~ 293 GeV, the value of the Hubble parameter is just H(M F ) = 1.2 x 1CT 4 eV, unable 
to excite even the lightest neutrino. However, the most interesting scales are much smaller. Let us 
first consider the nucleosynthesis epoch when the radiation dominates over the matter, and derive 
the restriction on the parameter v. In the RD regime, the solution for the density Q can be 
rewritten in terms of the temperature 



PR(T) = ^ (r-TY - [r^ - r 2 j Q% p° , (14) 

where r = T/Tq, Tq ~ 2.75 K = 2.37 x 10~ 4 eV is the present CMB temperature, g* = 2 for 
photons and g* = 3.36 if we take the neutrinos into account, but this difference has no importance 
for the present considerations. It is easy to see that the size of the parameter v gets restricted, 
because for v > 1 the density of radiation (in the flat case) would be the same or even below the 
one at the present universe. On the other hand near the nucleosynthesis time we have T 3> To, 
and so the corresponding value of the CC is 



Ar(T) ~ — ^— (r' u T) 4 



l-u 30 v ' l-2v 



. r 4(l-") _ r 2 



n% pi . (15) 



Hence, in order not to be ruled out by the nucleosynthesis, the ratio of the CC and the energy 
density at that time has to satisfy 

\^r/Pr\ ^ \u/(l-u)\ ~ \v\ <1. (16) 

A nontrivial range could e.g. be < \i>\ < 0.1. Both signs of v are in principle allowed provided the 
absolute value satisfies the previous constraint. Let us notice that, in view of the definition (|13|) . 
the condition v <C 1 also means that M < Mp. Hence, the nucleosynthesis constraint coincides 
with our general will to remain in the framework of the effective approach. It is remarkable that 
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Figure 1: Relative deviation 5H(z;u), Eq. (|2()|). of the Hubble parameter versus the redshift for 
various values of the single free parameter v in flat space with n° M = 0.3 and Q° A = 0.7. 

the two constraints which come from very different considerations, lead to the very same restriction 
on the unique free parameter of the model. The canonical choice M 2 = M p , corresponds to 

1 9 
\v\ = i/ ee ~ 2.6 x 10~ 2 (17) 

12 7T 

which certainly satisfies i/<1. Since the value of the parameter v is small, all the effects of the 
renormalization group running of the CC will be approximately linear in v . 



Some phenomenological Implications. 

Let us now turn to the "recent" Universe characterized by the redshift interval < z < 2, and 
let us evaluate some cosmological parameters which can be, in principle, extracted from the future 
observations, say by the SNAP project j^. The first relevant exponent is the relative deviation 
5 a of the CC from the constant value Aq. Then, using our solution ()10l) . we obtain (expanding at 
first order in v) 

A(*)~A + ^ [(l + ^) 3 -l] • (18) 
What about the numerical effect? The relative correction to the CC at redshift z is given by 

h = Afe %- A ' -.>gf [p + »y-i]. as) 

Let us take flat space with Q^ M = 0.3, Q A = 0.7, and the value v = vq defined in (|17|> . For z = 1.5 
(reachable by SNAP .1]) we find 5a = 16.3%, and so one that should be perfectly measurable by 
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Figure 2: Relative deviation 8Vtj^{z; v), Eq.(|21jl. of the cosmological constant parameter versus the 
redshift for the same cases as in Fig.^ 

SNAP. For larger values of v the test can be much more efficient of course. For instance, values 
of order v = 0.1 are still perfectly tenable, in which case the previous correction would be as large 
as > 60%! In general, the strong cubic z-dependence in 5\(z\v) should manifest itself in the 
future CC observational experiments where the range z > 1 will be tested. It is important to 
emphasize that v is the unique arbitrary parameter of this model for a variable CC. Therefore, 
the experimental verification of the above formula must consist in: i) pinning down the sign and 
value of the parameter v; and ii) fitting that formula to the experimental data. 

Next we present the relative deviation of the Hubble parameter H(z, v) with respect to the 
conventional one H(z, v = 0). At first order in is, 

H(z; v) - H(z; 0) 



5H{z;u) = 



H(z;0) 



^ M i + nS,[(i + *)3-i] • (20) 

Equation l|20|) gives the leading quantum correction to the Hubble parameter (JJJ) when the renor- 
malization effects in © and (|10j) are taken into account. Notice that 5H(0; v) = 0, because for 
all v we have the same initial conditions. Then for z / we have e.g. <5-ff (1.5; vq) ~ —2% and 
5H(2;i/q) ~ —3%. For higher v we get sizeable effects: 5H(z; 0.1) ~ —8% and —10% for z = 1.5 
and z = 2 respectively. See Fig. ^ for the detailed numerical evaluation. 

Although the induced corrections on H are not very high, the relative deviation of the renor- 
malized cosmological constant parameter VL^(z; v) = 8irGA(z;i')/3H 2 (z;i i ) with respect to the 
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standard one turns out to be much more significant. We get at leading order in u, 

ao r*.,A - n ^ u ) -Aa(z;0) , on 

<^M = n^o) W 

n° M (1 + zf - 1 1 + 3 n° M (1 + zf ln(l + z) 



n% n% + n% (i + z )3 

Again 5^(0; z^) = 0, as it should. Moreover, the deviation has the two expected limits for 
the infinite past and future, viz. 5£l\(oo;v) = oo and 5Q\(— l]v) = 0. The numerical evaluation 
of (|21|) is given in Fig. 13 Notice that even for v as small as vq, (|17j) . there is a sizeable 20% 
increase of £l\ at redshift z = 1.5 - reachable by SNAP. For v = 2vq the increase at z = 1.5 is 
huge, ~ 40%. For v < the effects go in the opposite direction. For fixed v, the larger is z the 
larger are the quantum effects. Thus if v ~ fo an d some experiment in the future can reach the 
far z = 2 region with enough statistics, then the effects on £l\(z) and the Hubble parameter can 
be dramatic. If they are seen, this model can provide an explanation for them. At present Oa 
has been determined at roughly 10% from both supernovae and CMB measurements, and in the 
future SNAP will clinch Vt\ to within ±0.05 The previous numbers show that for z > 1 the 
cosmological quantum corrections can be measured already for a modest v > 10~ 2 . A complete 
numerical analysis of this kind of FLRW models, including both the flat and curved space cases, 
together with a detailed comparison with the present and future Type la supernovae data, will be 
presented elsewhere ^] ■ 

Conclusions 

To summarize, we have presented a semiclassical FLRW type of cosmological model based on 
a running CC with the scale [i = H. We have shown that if the decoupling quantum effects on A 
have the usual form as for the massive fields, then we can get a handle on problem II by giving an 
alternative answer to problem III. In fact, we can explain the variation of the CC at low energies 
without resorting to any scalar field. The CC is mainly driven, without fine tuning, by the "relic" 
quantum effects from the physics of the highest available scale (the Planck scale), and its value 
naturally lies in the acceptable range. It is remarkable that all relevant information about the 
unknown world of the high energy physics is accumulated into a single parameter v. Finally, we 
have shown that the next generation of supernovae experiments, like SNAP, should be sensitive 
to v within its allowed range. A non-vanishing value of v produces a cubic dependence of the CC 
on z at high redshift, which should be well measurable by that experiment, if it is really there. 
Therefore an efficient check of this alternative framework, which might be the effective behavior 
common to a large class of high energy theories, can in principle be performed in the near future. 
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